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ABSTRACT 

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. 
Alfven waves and magneto-sonic waves are particular classes of MHD waves. These 
wave modes are clearly different and have pure properties in uniform plasmas of in- 
finite extent only. Due to plasma non- uniformity MHD waves have mixed properties 
and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is 
a quantity unequivocally related to Alfven waves as compression is for magneto-sonic 
waves. Here, we investigate MHD waves superimposed on a one-dimensional non- 
uniform straight cylinder with constant magnetic field. For a piecewise constant den- 
sity profile we find that the fundamental radial modes of the non-axisymmetric waves 
have the same properties as surface Alfven waves at a true discontinuity in density. 
Contrary to the classic Alfven waves in a uniform plasma of infinite extent, vortic- 
ity is zero everywhere except at the cylinder boundary. If the discontinuity in den- 
sity is replaced with a continuous variation of density, vorticity is spread out over the 
whole interval with non-uniform density. The fundamental radial modes of the non- 
axisymmetric waves do not need compression to exist unlike the radial overtones. In 
thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves 
with phase velocities between the internal and the external Alfven velocities can be 
considered as surface Alfven waves. On the contrary, the radial overtones can be re- 
lated to fast-like magneto-sonic modes. 



Subject headings: Sun: oscillations — Sun: corona — Sun: atmosphere — Magneto- 
hydrodynamics (MHD) — Waves 
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1. Introduction 



Alfven waves were predicted 70 years ago by H. Alf ven in a short paper entitled "Existence 
of Electromagnetic-Hydrodynamic Waves" ( Alf yenl 1 1 942h . Alf yen's predictio n was initially met 
with disbelief and was accepted only years later (|Falfhammar & Desslerll2006|) . Alfven waves are 
ubiquitous in magnetized pla smas, in fusion plasm a physics, geophysics, astrophysics and solar 



physics (see recent review by iGekelman et al. 



2011). 



Alfven waves are a particular class of magnetohydrodynamic (MHD) waves. MHD waves 
have become a subject of intense research in solar physic s, largely because observations clearly 



show that they are ubiquito us in the solar atmosphere (e.g. iTomczyk et al .1120071 : be Pontieu et al 



20071: Mcintosh et alJ201 II) . A pioneering theory paper on MHD waves is the paper by lEdwin & Roberts 
(11983b on MHD waves on an axisymmetric magnetic cylinder of piecew ise constant de nsity and 
consta nt straight magnetic field (see also, e.g., Iwentzellll979bl : Ispruitll982h . The paper by lEdwin & Robertsl 
(|1983l) has paved the way for new research as observations brought to light new information on 
MHD waves that was u nknown in the 1980s. For MHD waves in the solar corona Figure 4 of 
Edwin & Robertsl (| 19831) is often used as reference frame. In this figure a variety of MHD waves 
are displayed. Fast and slow waves are present but Alfven waves are apparently absent from the 
diagram. The question then arises, where are the Alfven waves? 

Compression and vorticity are key quantities for characterizing MHD waves. In a uniform 
plasma of infinite extent MHD waves have non-zero compression and zero vorticity in case of 
magneto-sonic waves, or zero compression and non-zero vorticity in case of Alfven waves. The ba- 
sic characteristics of the classic Alfven wave are that its motions are vortical and the total pressure 
in the plasma remains constant during the passage of the wave. For an inhomogeneous medium, 
however, the total pressure, in general, couples with the dy namics of the motion, and the assump- 



tion of neglect of pressure perturbations becomes invalid (|Hasegawa & Uberoilll982|) . However, 
that does not mean that the concept of Alfven waves is obsolete. In general in an inhomoge- 
neous plasma, MHD waves have mixed properties which can be traced back to the properties 
of the classic slow, fast, and Alfven waves in a homogeneous plasma of infinite extent. The 
degree to which the classic properties are present in a given MHD wave depends on the back- 



properties or coupled waves can lead to damping and was discussed by, e.g., 


Chen & ] 


iasegawa 


(1974); 


Tirrv & Goossens 
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1996): Goossens & de Groof (2001): Goossens et al. 


(2002a 


b. 2011): 


De Groof & Goossensu 


2 




>); 


Terradas et al. (2008b): Callv & Goossens (2008 


): 


Pascoe et al. (I2OIO. 


2011): 


Callv & Andries 
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2010 


), among many others. 



An MHD wave on an axi-symmetric one-dimensional cylindrical plasma is characterized by 
two wavenumbers, i.e., the azimuthal wavenumber, m, and the axial wavenumber, k z . In addition, 
modes can have different nodes in the radial direction and this number of nodes can be used to 
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further classify the modes. The term fundamental radial mode used here refers to waves that have 
no nodes in the radial part of the eigenfunction. The main objective of this paper is to show that the 
fundamental radial modes of the non-axisymmetric MHD waves with phase velocities between the 
internal, v A ,i, and the external, v A . e , Alfven velocities are surface Alfven waves. These modes were 
originally called fast waves by Edwin & Roberts! (119831). and the adjecti v e 'fast' became widely 
used in papers that followed the nomenclature of Edwin & Robertsl (11983b . iGoossens et al.l (|2009|) 
investigated the forces that drive these waves and found that t he magnetic tension force always 
dominates the pressure force. In addition IGoossens et al.l (120091) showed that compression is small 
in the particular case of thin tubes. Hence these waves have not the typical properties of fast 
magneto-sonic waves and behave more as Alfven waves. This lead IGoossens et al.l d2009l) to call 
these waves A lfvenic. The ad jective Alfvenic was used in order to point out that pu re Alfven waves 



as dec ribed by lAlfvenl (|l942|) can only exist in a uniform plasma of infinite extent. IGoossens et al 



(1201 If) reconsidered these waves in their section on quasi-modes and decided to call them surface 
Alfven waves. 

In the present paper we continue the theoretical investigation of the nature of the waves. We 
are concerned with vorticity since it is a quantity that is unequivocally related to Alfven waves. 
Pure Alfven waves in a homogeneous medium are the sole waves that represent vorticity perturba- 
tions and any spatial distribution of vorticity can be represented by means of Alfven waves. Here 
we shall argue clearly that, in a piecewise homogeneous model, whether it is in the planar or cylin- 
drical case, the Alfven surface waves have zero vorticity everywhere except at the discontinuity 
where all vorticity is concentrated. They are hence appropriately termed 'Alfven' surface waves 
and complement the pure Alfven waves which are bound to the separate homogeneous regions. In 
view of their properties, the fundamental radial eigenmodes in a magnetic cylinder with mi^O and 
phase velocities in the r ange VA,i to VA, e can be considered as surface Alfven waves. As a matter 
of fact. IWentzell (| 1979bl) was more ambitious and called all MHD waves with phase velocities be- 
tween v A ,i and v A , e surface Alfven waves. However, since the radial overtones have an oscillatory 
spatial behavior in the internal plasm a we prefer to keep the term surface to the fundamental radial 
modes only. In addition we note that llonsonl (| 1978b called th e surface wave s Alfvenic probably to 
distinguish them from the pure Alfven waves as described by lAlfvenl (| 19421) . 



The current paper is relevant to the recent disc ussion in the solar physics community on the 



nature of the observed ubiquitous, transvers e waves (fTomczyk et al.ll2007l : iDe Pontieu et al.ll2007 



Okamoto et al.ll2007l : Irvlclntosh et al.ll201 ll) . The original authors of these pa pers all claimed the 
detect ion of the Alfven waves in the solar corona. This claim was challenged by I Van Doorsselaere et al 
(I2008h who argu ed that the waves s hould be interpreted as fast kink (m = 1) wav es following the 
nomenclature of Edwin & Robertsl (|1983f) . The results of IGoossens et al.l (120091) and the present 
paper show that these waves are actually surface Alfven waves. 
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This paper is organized as follows. In Section |2] we briefly review the properties of Alfven 
waves and slow and fast magneto-sonic waves in a uniform plasma of infinite extent with a con- 
stant magnetic field. Then we try to understand how the properties of the MHD waves known 
for uniform plasmas of infinite extent are modified when the plasma is no longer uniform and/or 
confined to a finite volume. First, in Section |3] this is done for the case of surface Alfven waves 
in a true discontinuity in the Alfven velocity. Later, in Section |4] we extend our investigation to 
the case of MHD waves in non-uniform magnetic cylinders and show that the properties of the 
radially fundamental non-axi symmetric transverse waves in cylinders are remarkably similar to 
those of surface Alfven waves in a true discontinuity. Finally, Section [5] contains our discussion 
and the relevant conclusions of this work. Those readers interested in our results but that do not 
wish to go in detail through the mathematical derivations will probably find useful the summary 
and discussion of Section [5J 



2. MHD waves in a uniform plasma of infinite extent 

The basic equations for the discussion of linear ideal MHD waves superimposed on a static 
plasma are 

p— = -Vp' + -(VxB')xB + -(VxB)xB', 

Ot jX 11 

dB' 

— = Vx(vxB), 
ot 

d P' 

— + v • Vp = -ypV ■ v, 
ot 

dp' 

■f+v-Vp = -pV-v, (1) 
ot 

where p, p, and B are the equilibrium plasma pressure, density, and magnetic field, respectively. 
We take p and p uniform and B straight and constant, namely B = Bl z . In addition, p' , p', v, 
and B' are the Eulerian perturbations of the plasma pressure, density, velocity, and magnetic field, 
respectively, y is the adiabatic index, and ii is the magnetic permittivity. Equations (H)) assume that 
there are no equilibrium flows. For incompressible motions we have to replace the 3th equation of 
Equations (PQ) with 

V • v = 0, (2) 

and treat the plasma pressure perturbation, p', as an unknown function in addition to the three 
unknown components of velocity. 

In what follows we shall find it convenient to use the Lagrangian displacement, so that for 
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a static plasma up to linear order 



v = 



dj_ 

at' 



(3) 



and the Eulerian perturbation of total pressure, P' = p' + p' m , where p' is the perturbation of the 
gas pressure and p' m = B • W/p. is the perturbation of the magnetic pressure. 

With the use of the Lagrangian displacement, £, one may integrate the induction equa- 
tion, the energy equation and the continuity equation at once and use these expressions to 
eliminate all other variables from the momentum equation so that it takes the form: 



d 2 e 



(4) 



Here F is the force oper ator as derived first by lBernstein et alj (119581) . (see also lFrieman & Rotenberg 
1960l : lGoedbloedlll983l) . 



We consider a uniform plasma of infinite extent. We study linear planar harmonic waves and 
put the wave variables proportional to the factor 



exp(z'k • r - icot) 



(5) 



where k = (k x ,k y ,k z \ is the wave vector, r = (x,y,z) is the position vector, and u> the freque ncy. 
Variables that catch the basic physics of the waves (see, e.g., Thompson! 1962 : Goossens 2003 ) are 
the component of the displacement parallel to the equilibrium magnetic field, the compression, 
Y, and the component of vorticity along the magnetic field, Z. Y and Z are defined as 



Y = k-f = -iV-£ 

Z = (kX£)-l; = -/(VX£) ; . 

With the use of these variables the equations that govern the waves are 

co 2 £ z -k z v 2 Y = 0, 
k 2 v 2 A k z t z + (oj 2 -k 2 (v 2 s +v 2 A ))Y = 0, 
(a) 2 -ur A )z = 0. 



(6) 



(7) 



Here k 2 = k^+k 2 +k 2 , v 2 = yp/p is the square of the speed of sound, v 2 A = B 2 /pp is the square of the 
Alfven velocity, and lo a = k,VA is the Alfven frequency. The system of Equations d7J) is decoupled 
into two subsystems. The first two equations involve the variables £ z and Y but not Z. They are 
the equations that define the (slow and fast) magneto- sonic waves. These waves have compression 
and a component of velocity parallel to the magnetic field, but no vorticity. The third equation only 
contains Z. It defines the Alfven waves. Alfven waves have vorticity but no compression and no 
component of velocity parallel to the magnetic field. 
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2.1. Alfven waves 



The dispersion relation of Alfven waves is obtained from Equations © with £, = 0, Y = 0, 
and Z ^ 0. Thus the dispersion relation of Alfven waves is 



An Alfven wave exists for any wave vector k but its frequency only depends on the component 
of k parallel to the equilibrium magnetic field. The frequency is degenerate with respect to the 
components k x and k y of the wave vectors in planes normal to the equilibrium magnetic field 
lines. Since frequency only depends on k z , the group velocity, v gr , is always directed along the 
equilibrium magnetic field and is equal to va, while the phase velocity, v p h, is by definition directed 
along the wave vector k, namely v P h = va cos(#)lk, where 6 is the angle that k makes with B, and 
l k and 1 B are the unit vectors in the direction of k and B, respectively. Alfven waves are anisotropic 
in the extreme. We recall that in our discussion so far no particular choice has been made for k and 
the wave vector is general. 



The dispersion relation for the slow and fast magneto-sonic waves is readily obtained by 
imposing that the subsystem in Equations © for the variables £ and Y has a non-trivial solution. 
The result is 



where the + sign corresponds to the fast/slow magneto-sonic waves, respectively, and co c is the cusp 
frequency. Subscripts "f" and "si" denote quantities related to fast and slow waves, respectively. 
Let ft denote the ratio of the plasma pressure to the magnetic pressure. Then in the approximation 
/? — > (v s = 0), the slow waves disappear from the scene, i.e., co s \ = 0. Fast waves still remain in 
the zero-/? approximation with frequency a>f = kv A . On the contrary, for v s — » oo fast waves are 
removed to infinite frequencies and a particular form of slow waves remain with co s \ = a>\. 



CO = U)a = k z V\. 



(8) 



2.2. Magneto-sonic waves 




(9) 



2.3. Forces and motions 



It is instructive to write the equation of motion as 



dt 2 



-V h /" + T h , 



d% dp' 

»w - ~k- (10) 

where T h is (the horizontal component of) the magnetic tension force and Vh and £ h = (£ x , g y , 0)' are 
the gradient operator and the displacement, respectively, in horizontal planes perpendicular to the 
constant vertical magnetic field. The component of the displacement parallel to the magnetic field, 
is solely driven by plasma pressure and unaffected by magnetic forces. With the temporal and 
spatial dependence specified by Equation ©, and k h = (k x , k y , 0)' Equations (flOl) can be written as 

-pu?£ h = -ik h P' + T h 

T h = -pofeb, 
pco 2 % z = ik z p'. (11) 

The force in horizontal planes due to total pressure, II, is 

Il = -ik h P' = -p(co 2 -co 2 A )t h . (12) 

Hence the ratio of the horizontal components of total pressure force to the corresponding compo- 
nents of the magnetic tension force is 

2 2 ? 
- CO - CO . CO 

A(co 2 ) = T A = — -1. (13) 

Equation (PT3T) is general and applies to Alfven waves with co 2 = co 2 A , fast waves with co 2 = to 2 and 
to slow waves with co 2 = co 2 v 

For the Alfven waves A(co 2 A ) = so that the only restoring force is the magnetic tension force. 
The displacements are incompressible and confined to horizontal planes since = 0. For the 
Alfven waves the displacement £ = £ h and is perpendicular to the horizontal wave vector k h so as 
to make compression zero and vorticity maximal. Alfven waves are highly anisotropic and totally 
insensitive to the value of the sound speed of the plasma. 

On the contrary, magneto-sonic waves involve both the total pressure force and the magnetic 
tension force. It is now appropriate to consider the relative importance of horizontal compres- 
sion k h • £ h and longitudinal compression k z g z , which can be readily derived from the first line of 
Equation ©: 

k h -£ h _co 2 -k 2 v 2 s co 2 1 
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Equation (PT4l) is again applicable to fast waves with co 2 = to 2 and to slow waves with co 2 = co 2 v 
For Alfven waves with co 2 = (x> 2 A , it is largely irrelevant as both the denominator and the numer- 
ator vanish and no compression is involved at all. The equation tells us that the coupling of the 
longitudinal compression to the perpendicular compression depends on the sound speed. For fast 
waves the ratio is positive while it is negative for slow waves. This means that compression, Y 
(the sum of the denominator and numerator of (Q3)), is maximized for fast waves and reduced for 
slow waves. If the sound speed is small, the longitudinal dynamics decouples and is irrelevant to 
the perpendicular dynamics. The slow waves vanish and there are no longitudinal motions in the 
fast waves. On the contrary if the sound speed is large, the above ratio approaches -1, illustrating 
the perfect coupling in an incompressible medium. The fast waves are banned to infinity while for 
the slow modes lo\ — > co 2 A so that the dynamics in the perpendicular direction is dominated by the 
magnetic tension force. 



3. MHD waves on a true discontinuity 

In the previous section we have studied the properties of Alfven waves and magneto-sonic 
waves in a uniform plasma of infinite extent. In that configuration, Alfven waves are strikingly 
different from fast waves. 

Let us now see how a deviation from the uniform plasma of infinite extent adds new com- 
plexity. First of all consider a magnetic field that is still straight and unidirectional throughout 
space and let all equilibrium variation depend only on one cartesian coordinate x which is directed 
perpendicular to the equilibrium magnetic field. Now the wave variables can be put proportional to 
exp (iky.y + ik z z - itotj, but the functional dependence on x must remain unspecified. The relevant 
equations become: 



(oj 2 - 


<> 2 \ d ^ 


-- K 2 P', 


p(co 2 


dP' 

dx ' 




= 


= ik y P', 




- 


-- ik z v 2 s P', 



(15) 
where 

co" - (yj + vjXco 2 - co 2 c )(k 2 + k 2 ) 
(if + v* ' 



is a function of x. 



-9- 



We may combine the first two equations in Equations (1151) into a second order ordinary dif- 
ferential equation for P', namely 



1 



dP' 



= K 2 P'. 



(17) 



From Equations (1151) we can also compute the components of vorticity, V x in the z-direction as 



(VXf}-l Z = ikyP' — 

ax 



K p(aj 2 -co 2 A ) / 



(18) 



We shall use these equations in the following subsections. 



Here we will first consider an equilibrium consisting of two uniform plasmas separated by a 
sharp discontinuity in the Alfven velocity. The discontinuity coincides with the plane x = x . For 
simplicity, we take a constant magnetic field along the z-direction. Thus, the discontinuity in the 
Alfven velocity is introduced by a discontinuity in density as, 



p(x) 



Pi, if, x < xq, 

Pe, if, X > Xq, 



(19) 



where both p x and p e are constant s and p\ ■£ p e . The MHD wa ves of this plasma configuration 
have been studied before (see, e.g., Wentzel 1979a : Roberts|[l981 ). The aim here is to point out the 
differences between classic Alfven waves and surface Alfven waves. 

The basic Equations © are valid at both sides of the interface and K 2 (whenever it is negative) 
in Equation ( fTTT ) must be indentified with -k 2 x used earlier. The only difference is hence in the 
boundary conditions. In particular, the solutions at both sides of the discontinuity need to be 
matched together by ensuring continuity of both total pressure P' and normal displacement % x . The 
interaction between the two media allows for solutions in that range where K 2 > 0. In that case 
solutions are given by: 

P'(x) = Ai, e exp (±K Ue (x - Xq)) . (20) 

So whenever K 2 > the modes correspond to surface waves decaying away from the surface and 
their existence crucially involves the interaction between the two media. In contrast, the classical 
slow and fast waves discussed earlier correspond to solutions in the domain where K 2 < 0. 



3.1. Classic Alfven waves 



However, let us first retrieve the classic Alfven waves by requiring P' = and £ z = every- 
where. Hence also p' = and since V • £ = also p' = 0. We find two different solutions under 
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these conditions. The first solution has frequencies to 2 = to 2 Ai and motions satisfying 

for, x<x , 
(V X £) • \ z = 0, for, x > xq. 

The second solution has frequencies co 2 = u) 2 Ae and motions satisfying 

(Vx£)-L = 0, for, x<x , 
(Vxf)T,# 0, for, x> x . 

Hence, classic Alfven waves with different Alfven frequencies co Ai and u) 2 Ae are confined to the 
half spaces x < x and x > x , respectively. 



3.2. Surface Alfven waves 

Let us now turn to the motions that do involve P' 0. We have argued above that the classical 
slow and fast waves are related to the domain where K 2 < while if K 2 > the modes are surface 
waves relying on the interaction between the two media. The roots of K 2 as a function of frequency 
are easily determined as they correspond to the fast and slow frequencies © with k x = 0. Addition- 
ally K 2 changes sign at the cusp frequency. This allows to easily indentify the frequency ranges 
for the (classical) slow and fast modes and of the surface modes. Now consider two important 
limiting cases. In the incompressible limit the slow frequencies collapse to the Alfven frequency 
and the fast frequencies are banned to infinity, hence only surface type solutions remain and in fact 
K 2 » (k 2 +k 2 ). S imilarly, in what i s called the limit of 'nearly perpendicular propagation' (k 2 » k 2 ) 



where K 2 « k 2 . IWentzell (I1979al) already remarked that the incompressibility condition is a 
good approximation when k 2 » k\. Furthermore, the latter approximation is highly relevant to 
thin cylindrical tubes considered further in this paper. In both cases a solution to the dispersion 



relation expressing the matching of g x and P' at x = xq is found as: 

Pi +Pe 



Pi co A .+p e oJ Ae 

CO = ; = co±. (23) 



The frequency lies in between the internal and external Alfven frequency and hence by Equa- 
tion (fT3l) A can't be very large at neither side of the discontinuity. In fact, if the magnetic field is 
constant, -1 < A < 1 so that the perpendicular dynamics is always dominated by the tension force. 
One might be tempted to classify the wave as either fast or slow, since there is compression in- 
volved. In that case one would need to take the ratio defined in Equation ([141) into account. In 
the incompressible limit we would then necessarily be dealing with a slow mode while in a cold 
plasma it would need to be interpreted as a fast mode. However, in the limit of nearly perpendicular 
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propagation the perpendicular dynamics is completely insensitive to the sound speed. The mode is 
dominated by the perpendicular dynamics and the sound speed just influences to what extent this 
couples to longitudinal motions as well. The insensitivity to the sound speed and the dominance 
of tension as a driving force is a strong argument to call these modes Alfvenic. 



Let us now focus on vorticity. From Equation (1181 ) we can evaluate the z-component of vor- 
ticity of the surface Alfven wave. To do so we note that p (co 2 - toj^j is a piecewise constant in our 
equilibrium and can be expressed as 



1 



A + Pe 



1 



PiPe 0)\ A -0)\. 



(l-2ff (*-*>)), 



p(co 2 -co 2 A ) 

where H(x - xo) is the Heaviside function defined as 

„, s f 1, if, X > X , 

H(x - x Q ) = < „ 

{ 0, it, X < x . 

Thus it is straightforward to write the expression for (V x £) • l z as 

vA+Pe 1 



(24) 



(25) 



(V x £) • l z = -likyF 



PiPe CO 



6(x - X ), 



A,e 



CO 



(26) 



A.i 



where 8{x - x ) is the delta function defined as 



S(x - x ) = 



1, if, x — Xo, 
0, if, x ± xo- 



(27) 



Equation (1261) shows that the surface Alfven wave is different from the classic Alfven wave. For the 
classic Alfven wave (V x £) • 1- is different from zero everywhere in the appropriate half space x < 
xo or x > x . The surface Alfven wave has (V x g) ■ 1, = everywhere except at the discontinuity 

X = Xo- 

When the true discontinuity in va is replaced with a continuous variation from VA,i to VA, e 
then the interval [o>A,i, <^A,e] is filled with the continuous spectrum of shear Alfven waves. Each 
magnetic surface, i.e., a surface of constant v A , now oscillates at its local Alfven frequency. The 
frequency of the surface Alfven wave is in the Alfven continuous spectrum and gets damped by 



^-^^ . M1 „ „ , „ — r^^-^*., „„, » ^ . 

frequency equals the frequency of the surface wave, (see, e.g., Tataronis & Grossmann 


1973; 


Grossmann & Tataronis 1973: Ionsonl 1978: Haseaawa & Ubero 


i 1982 


: Goedbloedlll983 


Hollwea 


1987b a: Hollweg & Yang 1988: Ruderman & Goossens 1993. 


1996|: I 


luderman et al. 


1995 


). Vor- 



ticity is no longer confined to a single surface but is spread out over the whole region of non- 
uniformity as indicated by Equation (TTSl . 



4. MHD waves in a magnetic cylinder 



In the present section we are interested in MHD waves superimposed on a ID straight cylinder. 
First we consider a piece- wise constant density as in Edwin & Roberts! (| 1 9 8 3b . Later we replace the 
jump in density with a continuous v ariation of dens i ty. Th e properties of kink MHD waves in ID 
straight cylinders were discussed by lGoossens et al.l(|2009|) . As for the case of a true discontinuity, 
here we give special attention to vorticity. We shall show that the surface Alfven waves in a true 
discontinuity and the radially fundamental MHD waves with phase velocities between VA,i and VA, e 
in a cylinder have strikingly similar properties. 



4.1. General theory 



The equilibrium configuration is a straight cylinder. We use cylindrical coordinates, namely 
r, (p, and z, for the radial, azimuthal, and longitudinal coordinates, respectively. The equilibrium 
quantities are functions of r only. Since the background model is independent of the spatial coor- 
dinates <p and z, and of time, t, the perturbed quantities are put proportional to 



exp (i(mcp + k z z - cut)) , 



(28) 



where m is the azimuthal wave number, k z is the longitudinal wave number, and co is the frequency 
as before. Since the background is variable in the radial direction there is not really a radial wave 
number. Alternatively, we can use the number of nodes in the radial part of the eigenfunctions 
to distinguish between radial fundamental and overtone modes. The equations for linear MHD 
motions superimposed on a I D cylindrical equilibrium mo del can be found in, e.g., lAppert et al. 
(119741) : ISakurai et al.l (|199ll) : iGoossens et al.l (| 19921 Il995|) . For a straight and constant magnetic 
field the equations for linear displacements on a ID cylindrical equilibrium model take the follow- 
ing simplified form 



D 



d(r£) 

dr 
dP' 



p(oj 2 - oo 2 A )£ v 



p(oj 2 - 



= -C 7 rP', 



p(oj 2 - oj\)^ r , 



an , 

—P', 
r 



ik 



^ 2 + v 2 A 



P\ 



(29) 



where all quantities have the same meaning as in previous sections. We recall that now the equi- 
librium quantities are functions of r. The coefficient functions D and are 
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D = p(v 2 s +v 2 A )(to 2 -aj 2 A )(to 2 -to 2 c ), 



C 2 = aj 4 -(v 2 + v 2 A )((o 2 -to 2 c )(^-+k 2 ). 



(30) 



Since lAppert et al.l (|1974|) it is well known that the system formed by Equations ( 1291) has 
mobile regular singularities at the positions r = r A and r = r c where to = to A (r A ) and to = to c (r c ). 
This leads to the definition of Alfven and slow (or cusp) continuous parts, namely 



[min coa> max &>aL [min to c , max to c ] , 



(3D 



with singular eigensolut ions (see, e.g.. lGoedbloedlll983HSakurai et al.lll99lHGoossens et al.lll995 
Tirry & Goossenslll996|) . For a straight field the (p- and z-directions are the directions of constant 
Alfven velocity perpendicular and parallel to the magnetic field lines, respectively. The r-direction 
is normal to the surfaces of constant Alfven velocity. Hence for a straight field ^ = £ x is the 
characteristic quantity for the Alfven waves and £. = £n for the slow waves, where || and _L de- 
note the directions parallel and perpendicular to the equilibrium magnetic field, respectively. £ r 
characterizes the fast magneto-sonic waves. 

For an equilibrium with a straight magnetic field, the Eulerian perturbation of total pressure P' 
is the quantity that produces waves with mixed Alfven and magneto-sonic properties. The coupling 
function, Ca, is 

mB 

(32) 



mB , 

C A = —P'. 

r 



The coupling function C A was introduced by 



Sakurai et al 



(1991) andlGoossens et al.ldl992h. and 



its role was discussed by, e.g., iGoossensI (|2008|) and lGoossens et al.1 (|201ll) . 

The two first order differential equations in Equations (|29l can be rewritten as a second order 
ordinary differential equation for P' as 



p(oj 2 - co 2 A ) d 
r dr 



dP' 



p(to 2 - co.) dr 



nr 



p', 



where T(to 2 ) is defined as 



T(co 2 ) 



(to 2 - k 2 v 2 ,)(co 2 - to 2 A ) 
(vi + vDitoi-col) ' 



(33) 



(34) 



Note that Equation ([111) is applicable here also with £ h = % r l r + g v l v and k h = -t-Ajz + 
Equations (fT2)) and (TT3l) remain valid so that we have a simple expression for A to decide on the 
relative contributions of pressure and tension forces. 

Recall that the component of vorticity parallel to the equilibrium magnetic field, (V x if) ■ l z , 
is unequivocally related to Alfven waves. Alfven waves have (V x if) • l z + everywhere in an 
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infinite and uniform plasma, while surface Alfven waves in a true discontinuity have (V x£) • 1, ^ 
at the discontinuity only. Hence for comparison with these previous cases and for later use it is 
instructive to derive an equation for this quantity: 

(VX0-1, = i-P'U 2 l ). (35) 
r dr \p(u> 2 - coj) J 



Note that this equation is formally equivalent to Equation (PT8I ) derived for waves on a true discon- 
tinuity if in Equation (PT8"1) k y is replaced by m/r and the derivative in x is replaced by a derivative 
in r. 



4.2. Piecewise uniform plasma 



Here we consider the case studied by Edwin & Robertsl (|1983l) . i.e., the case of a piecewise 
constant density profile. The general situation in which the density varies continuously is consid- 
ered later. We assume the following density profile, 



P(r) 



pi, if, r < R, 
p e , if, r > R, 



(36) 



where both p y and p e are constants and R denotes the radius of the cylinder. Subscripts i and e refer 
to internal and external plasmas, respectively. We focus on the case p x > p e . Thus, there is a jump 
in density, namely p\ - p e , at the cylinder boundary. 

We rewrite Equations (1331 and (1351 when the density p and the local Alfven frequency a> A are 
both constants, namely 



1 d 

r dr 

p( W 2 -^)(Vx^).l : = 





r 2 i 
m 


ffl- 





P', 



(37) 
(38) 



Equation (T37T) is Equati on (5) of Edwin & Robertsl (|1983|) . Equation (T5HT ) is Equation (3b) of 
Edwin & Robertsl (|1983h when the dependency of Equation (T28l) for the perturbed variables is 
used. 

Notice that Equation (1331) and, a fortiori, Equation (137T > are very reminiscent of Equation (fTTT ) in the 
Cartesian geometry. In the incompressible limit T —* -kj and the analogy is clear if one identifies 
m/r with k y as done earlier. But more importantly, if one considers 'thin tubes', the right hand 
side coefficient is dominated by (m/r) 2 (in the domain of interest, i.e. r ~ R), which is analogous 
to the limit of 'nearly perpendicular propagation' considered in the Cartesian case. In particular 
we conclude again that Equation (1371) thus becomes ignorant of and insensible to the value of the 
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sound speed. Furthermore, the solutions decay away from the discontinuity surface, indicative of 
surface wave behaviour. Nodes are only found in the solutions at radial distances r » R, far away 
from the domain of interest. 



4.2.1. Classic Alfven waves 



Edwin & Robertsl (|1983l) concentrated on obtaining solutions to their Equation (5), our Equa- 
tion (1371) . for a piecewise uniform cylindrical plasma and did not consider their Equation (3b), our 
Equation (1381) . further. Let us now focus on the solutions to Equation (1381) and investigate what 
has happened to the classic Alfven waves when we move from an infinite uniform plasma to a 
cylindrical plasma with a piecewise constant density. 

First we consider the extreme case that p; = p e so that the plasma is uniform. In that case co 2 A 
is a constant and we have a solution of the system of Equations (1291 for any m, namely 

co 2 = co 2 A (39) 

with 



^ * o, &*o, £ = o, 

P' = 0, ^ = 0, V-f = 0. (40) 
dr 

The only constraint is that the waves described by Equations (I3"yi)-(l4"0"l) have to satisfy is V • g = 0. 
This can be done in many ways. The only restoring force is the magnetic tension force T = 
-pco 2 A {i; r \ r + ijtpltp) = -ptojj;. Thus, in an infinite and uniform cylindrical plasma pure Alfven 
waves are independent of the azimuthal number number m. 

We turn to p x p e . Again the solutions must satisfy the constraints given in Equations (T39l 
and (1401) . Since co 2 = co 2 Ai for r < R and co 2 = co 2 Ae for r > R we have Alfven waves that live 
in the interior and in the exterior, respectively, of the flux tube. For co 2 = co 2 A . the components g r 
and ^ are different from zero for r < R and are identically zero for r > R. In addition ^ r = at 
the boundary r = R. In the particular case m = 0, <; r = everywhere. Conversely for co 2 = co 2 A 
the components £ r and % v are different from zero for r > R and are identically zero for r < R. 
At the boundary r = R again f r = 0. Note that the frequencies are independent of the azimuthal 
wavenumber m. When we replace the piecewise constant density profile by a fully non-uniform 
density profile, the modes with m = are the only ones that survive as purely incompressible 
modes. 
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4.2.2. Surface Alfven waves 

Here we investigate the solutions to Equation (1371) with m ± 0. We vary p, - p e and see what 
happens to the part of the spectrum with phase velocities between vaj and v A , e ]. We start from 
the dispersion curves of solutions for p; - p e = l-5p e , because this is analogous to Figure 4 in 



Edwin & Robertsl (119831) . We keep p- x constant and decrease p ; - p e from 1 .5p e to and follow the 
evolution of eigenmodes on the original dispersion curve. This evolution is illustrated in Figure Q] 
for m = 1 modes, i.e., kink modes. 

In the top graph of Figure \T\ we show the dispersion diagram of compressible modes with 
the phase speed on the vertical axis and the normalized wave number on the horizontal axis, for 
different density contrasts in various colors. The associated external Alfven speed v A , e (dashed line) 
and the density contrast are also indicated. In the solution, the internal Alfven speed has always 
been normalized to v A ,i = 1. In the bottom graph of Figured] we display the equivalent dispersion 
diagram but for incompressible modes. We can compare the top and bottom graphs of Figure Q] 
to assess the differences between compressible and incompressible waves. On purpose we have 
plotted the bottom graph on the same scale as the top graph in order to make the differences as clear 
as possible. The more striking difference is that the upper part of the bottom (incompressible) graph 
is empty. The dispersion curves in the top right belong to radial overtones of the kink modes. For 
radial overtones the total pressure perturbation P' has an additional node in the internal region. We 
notice the absence of the radial overtones in the incompressible dispersion diagram. The dispersion 
curves in the bottom of both graphs in Figure [Dbelong to the fundamental radial kink mode. The 
fundamental radial modes survive in the incompressible limit while radial harmonics are absent. 
Radial harmonics need compression to exist. In contrast the fundamental radial modes do not need 
compression to exist. Compression is a typical characteristic of magneto-sonic waves. Hence the 
fundamental radial modes do not have the typical properties of fast magneto-sonic modes. Instead, 
they behave like surface Alfven waves. 

The behavior of the fundamental radial modes in the thin tube limit (k z R <sc 1) is the same in 
both compressible and incompressible cases, i.e., their phase velocity tends to the kink velocity, 
Vk, namely 

/PiV Aii +PeV A>e 

v k = \ : -• (41) 

\ Pi+Pe 

The behavior of compressible and incompressible waves is slightly different for large k z R, i.e., 
the phase velocity of the compressible waves tends to v A ,i while that of incompressible waves 
remains between v A ,i and the kink velocity, v^. We display in Figure |2] the dispersion diagram of 
the fundamental radial modes with different values of m and a fixed value of the density contrast. 
In the thin tube limit, i.e., k z R «; 1, the results for the different values of m overlap and their 
phase velocity tends to the kink velocity, v k . The differences between the modes with different 
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Fig. 1. — Top graph: Dispersion diagram a>/k z vs. k z R for compressible kink (m = 1) modes with 
different density contrasts plotted with different colors. The respective external Alfven velocity is 
displayed with a horizontal dashed line, also showing the density contrast. Bottom graph: Same as 
the top graph but for incompressible waves. 
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values of m grow when we take larger values of k z R. For k z R <sc 1 the frequencies are to a good 
approximation independent of k z R and independent of m. The behavior of the fundamental radial 
modes is reminiscent of Alfven waves in uniform plasmas of infinite extent and of surface Alfven 
waves in a true discontinuity. 

In Figure [3] the evolution of the phase speed with varying density contrast is shown for two 
different wave numbers, namely k z R = 0.3 and k z R = 2.3 (indicated as vertical dashed lines in the 
top graph of Fig. [Q. Several observations can be made from the results of Figure |3] It is clear 
that the dispersion curve for the fundamental radial kink mode is always between the external and 
internal Alfven speed. When these two values approach each other, i.e., pi -p e — > 0, the dispersion 
curve eventually collapses to the Alfven velocity, a>/k z = va,i- In that sense the fundamental radial 
mode is the descendant of the Alfven mode of the uniform case. In fact, we have calculated that 
all curves that start at the kink speed (Equation (I4TT )) collapse to the internal Alfven speed as the 
density difference p\ - p e goes to 0. The evolution of the radial overtones is entirely different from 
that of the fundamental radial mode. As the density contrast decreases, the radial overtones are 
less well confined to the magnetic cylinder. When the density contrast reaches a critical point, 
confinement is completely breached and energy leaks away in the form of MHD radiation. Then, 
radial overtones be come leaky modes with complex frequencies because of damping d ue to MHD 
radiation (see, e.g., lwilsonfl98 ll : Ispruitlll982[ ICallvlfl986T: iGoossens & Hollwegll 19931) . When the 
density contrast decreases, the phase speeds of these radial overtones do not collapse to the Alfven 
velocity. These modes are not related to the Alfven mode of the uniform case, but rather to fast 
modes. As such, they become leaky when the density contrast is too low. This same argument 
explains the peculiar behavior of the fundamental radial mode with m = 0, i.e., the sausage mode 
(not displayed in Figured). 



4.2.3. Vorticity 



Here we focus on the role of vorticity. As explained before, vorticity is a typical characteristic 
of Alfven waves. Let us now consider the vorticity of the solutions of Equation (T37T ). 



The reader of the original paper of Edwin & Robertsl (|1983l) might have the impression that 
vorticity for the solutions of a piecewise constant equilibrium vanishes. However, visual inspection 
of the radial variation of ^ shows that it is dis continuous at r = R w ith opposite values to the 
left and right of th e boundary (see Figure lb of IGoossens et al .1120091 corresponding to the m = 
1 mode). Actually iTerradas et al.l (|2008a|) studied the Kelvin-Helmholtz instability triggered by 
the velocity shear in ^ at the boundary. So there is vorticity present in this configuration. For 
mathematical simplicity let us adopt the thin tube (TT) approximation, i.e., k z R <sc 1, so that the 
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Fig. 3. — Evolution of the phase speed as a function of density contrast, for k z R = 0.3 (top graph) 
and k z R = 2.3 (bottom graph). The external Alfven speed is indicated with a dashed line, while the 
internal Alfven speed is set to 1. These two graphs are vertical cuts in the top diagram of Figure [T] 
along the indicated vertical long dashed lines. In the bottom graph n = means the fundamental 
radial mode and n = 1 means the first radial overtone. 
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dispersion relation for modes with m + is (see details in, e.g. jGoossens et all 19921 120091) 



Pi co 2 Al +p e co 2 Ae 

0) = ; = (O k . (42) 

Pi +Pe 

Note that this is exactly the same dispersion relation as for incompressible surface Alfven waves 
at a true discontinuity (Equation (1231). Equation (1421 is independent of m. 



From Equation (1421 it follows that 

1 A + Pe 1 



P(U) 2 -0) 2 A ) Pipe OJ Ae -CO Ai 



[l-2H(r-R)]. (43) 



Here H is the Heaviside step function as defined in Equation (1231) . Since P' is continuous at r = R 
it follows from Equations (1351) and (l43l that 

(V X {) ■ l z = -2i^P' Pi+Pe 2 1 2 8{r - R), (44) 

* PiPe < e -<i 

where 5 is again the Dirac delta function. Thus, the solutions to Equation ( 1371 ) for the fundamental 
radial mode of MHD waves with frequencies to e [coas, coA,e] have vorticity but it is concentrated 
as a delta function at the boundary. This is exactly the same behavior obtained for surface Alfven 
waves at a true discontinuity (Sect. I3.21 ). This is a rather pathological situation that finds its origin 
in the fact that the equilibrium has been forced to be piecewise uniform. The singularities of the 
continuous spectrum are all concentrated in the point r = R. Note that this is true for all modes 
with m ± 0. This result is a strong argument in favor of a classification of the fundamental radial 
modes as surface Alfven waves instead of fast body modes. 



4.3. Continuous density variation 

In the previous Subsections we have studied the properties of the waves for a piecewise con- 
stant density profile. Here we replace the discontinuity in density with a continuous variation 
of density in an intermediate layer of thickness /. Thus, density is non-uniform in the interval 
]R - 1/2, R + l/2[. Since co A (r) is non-constant in the interval ]R - 1/2, R + l/2[, it follows from 
Equation (1351) that (V x • \ z 4- in that interval. This means that vorticity is spread out over the 
interval with non-uniform density. 

At this point it is in structive to note that when we replace the piecewise constant density of 



Edwin & Robertsl (|1983|) by a continuous variation of density all wave modes with phase velocities 
between v A ,i and v A , e are in the Alfven continuous spectrum. As a consequence, the waves with 
m ± undergo resonant damping. The fundamental conservation law for resonant Alfven waves 
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was obtained by ISakurai et al.l (|199l|) in ideal MHD and by iGoossens et al.l (| 19951) in dissipative 
MHD for the driven problem, and by iTirry & Goossensl (|1996|) for the eigenvalue problem. For a 
straight magnetic field the conserved quantity is the total pressure perturbation, P' . 

We denote by r A the position of the Alfven resonant point where co = a> A (r A ) and use a Taylor 
expansion of co 2 - co\ in the vicinity of r A , namely 



CO 



CO, 



sAf 



(45) 



where A A = ^ (co 2 - co 2 ^j and s = r - r A . It follows from the third equation in Equations (|29l that 
in ideal MHD ^ diverges as l/s near the resonant point. Equation (|33T > implies that the singular 
behavior of (V x • l z is l/s 2 and hence stronger than that of Hence vorticity is different 
from zero everywhere in the non-uniform plasma when dco 2 A /dr ± 0, but it is by far largest at the 
resonant position r = r A . 

Now we can compare the behavior of vorticity in the non-uniform case with that in the par- 
ticular case of a piecewise constant density profile (Equation (1441)). The delta function behavior at 
r = R obtained for the piecewise constant profile is a very pathological and peculiar situation. In 
the non-uniform case vorticity is spread everywhere in the region of non-uniform density. 

Conversely for a wave with its frequency in the slow continuum, £ z and V • £ are singular at 
the slow resonant position r c where co = co c (r c ). Both quantities diverge as l/s with s now defined 
as s = r - r c . For coronal conditions the Alfven continuum and the slow continuum do not overlap. 
Hence when we study MHD waves with frequencies in the Alfven continuum we do not have to 
worry about the s low continuum. The situation is different in, e.g., thin threads of prominences (see 
Soler et al.l 120091) . and photopheric flux tubes in which the frequency of the radially fundamental 
modes with m ^ is in both the Alfven and slow continua. 

Next we compute the eigenfrequencies and the perturbations of the fundamental radial mode 
numerically. In order to have non- singular eigensolutions we need to remove the singularity by 
including dissipative effects. For that reason we compute eigenmodes of non-uniform equilibrium 
states in resistive MHD. We add the term ?7V 2 B' to the right-hand side of the linearized induction 
equation (the second equation of Equations (OQ)), where rj is the coefficient of magnetic diffusion 
or resistivity. For simplicity we take 77 as a constant. We define the magnetic Reynolds number as 
= v Ki R/r]. 



Our numerical procedure uses the PDE2D code (|Sewelll |2005|) based on finite elements to 
solve the eigenvalue problem defined by Equations (Q3 in our equilibrium. The numerical integra- 
tion of Equations (OQ) is performed from the cylinder axis, r = 0, to the finite edge of the numerical 
domain, r = r max , which is located far enough to obtain a good convergence of the solution and to 
avoid numerical errors. This means that we take r max » R. We use a nonuniform grid with a large 
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Fig. 4. — Top: Absolute value of the z-component of vorticity (in dimensionless units) in the non- 
uniform layer for the fundamental radial m = 1 mode with k z R = 0.1, Pi/p e = 2, and l/R = 0.2. 
The solid line is the resistive result with R m = 10 7 . The dashed line shows the ideal spread of 
vorticity due to non-uniformity computed from Equation (l3~5l) using the numerically obtained P' . 
Bottom: Same as the top panel for the resistive result but with different values of l/R. 
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density of grid points within the nonuniform interval ]R - 1/2, R + l/2[. The nonuniform grid also 
allows us to correctly describe the small spatial scales of the eigenfunctions within the nonuniform 
region due to the Alfven resonance. The PDE2D code uses a collocation method and the general- 
ized matrix eigenvalue problem is solved using the shifted inverse power method. The output of 
the program is the closest complex eigenvalue to an initial provided guess and the corresponding 
perturbations. 

First we use the components of the displacement, obtained numerically, to compute (V x£) • l z 
in the non-uniform region for l/R = 0.2 (see Figure@k). The remaining parameters are given in the 
caption of the Figure. Since (V x • l z is a complex quantity, we plot its absolute value. We use 
dimensionless units so that the maximum of |(V x • 1,| has been set to unity. Vorticity is present 
in the whole non-uniform region and is maximal at the Alfven resonance position, r A « R. We 
overplot in Figure Bit the ideal |(V x ■ l z \ computed from Equation (|35l ) using the obtained P' 
from the numerical code (see the dashed line). From either curve we conclude that non- uniformity 
spreads out vorticity but the resonant behavior is so strongly present that the values of (V x • 1, 
close to the ideal singularity totally overpower the values away from that position. 

Equation (|35l) gives us the ideal behavior of vorticity and doesn't include the effect of diffu- 
sion. Diffusion removes the singular 1 / s 2 behavior of (V x • 1, found in ideal MHD and limit 
(V x • 1- to finite values. The remnant of the ideal 1 / s 2 behavior is, however, still clearly present. 
In fact, the effect of diffusi on is important in a dissipative layer of width <5 A around the resonance 



position given by (see, e.g. JSakurai et al.lll99l|) 



For the particular case of Figure Ek, the dissipative layer extends approximately in the interval 
0.95 < r/R < 1.05. The resistive result of Figure corresponds to R m = 10 7 , while the actual 
Reynolds number in the corona is believed to be around R m = 10 12 . Using the actual value of R m 
is unpractical from the computational point of view as it requires taking an enormous number of 
grid points in the numerical domain. We therefore use a smaller value of R m , and consequently 
a smaller number of grid points, but the qualitative effect of magnetic diffusion remains correctly 
described. Thus, the width of the dissipative layer around the resonance position would be ex- 
tremely thin if the actual Reynolds number would be used. We have made sure that the Reynolds 
numbers used in the computations are in the so-called plateau where the damping rate by resonant 
absorption is independent of R,„ and so the wave beh avior is not dominated by diffusion (see, e.g., 



Poedts & Kernerll 19911 : IVan Doorsselaere et al.l l2004) 



Now we vary the thickness of the non-uniform region, /. Figure |U> displays vorticity in the 
nonuniform layer for different values of / /R. For comparison purposes we have set max |(V x • l z \ 
1 in all cases. As before, vorticity is larger near the resonance position, but vorticity spreads along 
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the whole region of non-uniform density. In a fully non-uniform equilibrium, vorticity would 
spread out over the whole domain. 



5. Discussion 

In this paper we have studied the different properties of linear Alfven waves and magneto- 
sonic waves in uniform and non-uniform plasmas. First in a uniform plasma of infinite extent, we 
have reiterated that Alfven waves are driven solely by the magnetic tension force and that they are 
the only waves that propagate vorticity. The displacements are vortical and incompressible. On the 
contrary, magneto-sonic waves are driven by both the total pressure force and the magnetic tension 
force. The displacements are compressible and have no vorticity. 

Then we have moved to non-uniform plasmas and have investigated how MHD waves are 
modified by non-uniformity. For the case of a true discontinuity in the Alfven velocity we find 
that the incompressible surface Alfven waves have vorticity equal to zero everywhere except at 
the discontinuity, where all vorticity is concentrated. The behavior of the surface Alfven waves is 
clearly different from that of the classic Alfven waves in a uniform plasma of infinite extent, which 
have vorticity different from zero everywhere. 

Subsequently we have considered the case of MHD waves superimposed on a ID non-uniform 
straight cylinder with const ant magnetic field . For the particular case of a piecewise constant 



density profile as studied by Edwin & Robertsl (11983b . we find that the fundamental radial modes 



of the non-axi symmetric (m 4- 0) waves with phase velocity between VA,i and VA, e have properties 
remarkably similar to those of surface Alfven waves in a true discontinuity. In this pathological 
situation vorticity is present as a delta function at the cylinder boundary. When the discontinuity 
in density is replaced with a continuous variation of density, vorticity is spread out over the whole 
interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves 
do not need compression to exist unlike the radial overtones. 

With th ese insights we may now also inter pret the physics behind the computational results 



presented in I Van Doors selaere & Poedtsl (|2007|) . In that article the evolution of the m = 1 kink 
mode frequency was followed while the thickness of the inhomogeneous layer around the flux tube 
was increased (l/R /'"). It was found that the kink mode frequency joined the Alfven continuum 
when l/R passed a critical threshold. Indeed, understanding these waves as surface Alfven waves, 
now lets us conclude that the frequencies re turn to the Alfv en continuum as the surface 'goes 



away'. No surface, no surface mode. See also ISedlacekl (| 197 1|) for a similar interpretation 



We would like to stress that the importance of the labels Alfven or fast is not in the names 
themselves but in the properties that are intrinsically associated with these names. Due to plasma 
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non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or 
pure magneto-sonic waves. However, there are basic characteristics that remain strongly related 
to the wave type. Our results show that in ID magnetic cylinders the fundamental radial modes 
of the waves with m and phase velocities between v A ,i and v A , e have not the typical properties 
expected for fast magneto-sonic waves. Instead, their properties resemble very much those of 
surface Alfven waves in a tru e discont i nuity. F or this reason we call these waves surface Alfven 
waves as was already done by iWentzell (|1979bl) . 



Here we go back to t he discussion on the nature of the ubiquitou s , transverse waves as ob- 
served in the solar corona (ITomczvk et al.ll2007t iDe Pontieu et al.ll2007l : iMcIntosh et al.ll201 lb . In 
view of the results of the paper by IGoossens et al. ( 2009 ) and the present paper, the controversy 
about the interpretation of the observed waves is lifted. The fundamental radial modes of kink 
(m = 1) waves with phase velocity between the internal and external Alfven v elocities can be 



consi dered as surface Alfven waves (or Alfvenic waves in the nomenclature of iGoossens et al 



20091) . The two interpretations refer to the same physical pheno menon of a wave dominated by 
tension forces. The controversy was also partly due to the claim by IVan Doorsselaere et al.l (|2008|) 
that Alfven waves have to be torsional, i.e., axi-symmetric with azimuthal wavenumber m = 0. 
Axi-symmetric MHD waves do not displace the axis of the magnetic cylinder and the cylinder 
as a whole. The view that Alfven waves need to be axi-symmetric is too narrow. Anyway, in a 
non-uniform plasma there is a continuous spectrum of Alfven waves with frequencies independent 
of the azimuthal wavenumber, m. 



Note that the observations of ITomczvk et al.l (|2007l) are not the first observations of Alfvenic 
waves, but that they are the first to observe the ubiquity of these waves. Surface Alfven waves as de- 
scribed in the present paper have been observed on earlier accou nts although the author s at that time 
did not realize that they had indeed observed Alfvenic waves. IGoossens et al.1 (|2009l) pointed out 
that accepting resonant absorption as damping mechanism of the transverse MHD waves observed 
with the Transition-Region And Corona l Explorer (TRACE) im plied that these waves are surface 
Alfven waves or, in the nomenclature of IGoossens et al.l (|2009|) . Alfvenic waves. The fundamental 
radial modes of kink waves with phase velocity between the internal and external Alfven veloci- 
ties are surfa ce Alfven wave s . Hen c e the TRACE observation s of t ransverse MHD waves st arting 
in 1999 with lSchrijver et al.l (|1999l) . Aschwanden et al.l (|1999|) . and lNakariakov et al.l (| 19991) were 
observations of surface Alfven waves. 

We finally note that although the observed waves can be interpreted as a type of Alfven waves 
it is crucial to consider an adequate description of the coronal loop structures in which these waves 
propagate. Hence the use of expressions for Alfven waves in uniform plasmas of infinite extent 
may be inadequate for the study of waves propagating in the solar corona. A description in terms of 
surface Alfven kink (m = 1) waves in cylinders may offer a more detailed model for the seismology 
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and calculation of the energy budget of the observed waves. 
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